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ps^ , Summary. We give a geometrical characterization of A-prolongations of vector 

.^ ' fields, and hence of A-symnietrics of ODEs. This allows an extension to the case 

f^ , of PDEs and systems of PDEs; in this context the central object is a horizontal 

one-form /x, and we speak of //-prolongations of vector fields and //-symmetries 
O , of PDEs. We show that these are as good as standard symmetries in providing 

I ' symmetry reduction of PDEs and systems, and explicit invariant solutions. 

Introduction 



^ 



X 



The study of differential equations was the main motivation leading S. Lie to 
create what is now known as the theory of Lie groups; symmetry methods for 
; ^ ■ differential equations have received an ever increasing attention in the last fifteen 

C^ \ years, and by now there is an extremely vast literature devoted to these and/or 

their applications. 

It was recently pointed out by Muriel and Romero [8] that, beside standard 
symmetries, another class of transformations is equally useful in providing sym- 
metry reduction for scalar ordinary differential equations (ODEs); these were 
christened C°° symmetries^ or even X-sym,m,etries, as they depend on a smooth 
scalar function A (see also [5, 9] for applications of A-symmetries) . Soon after- 
wards, Pucci and Saccomandi identified the most general class of transforma- 
tions sharing the "useful" properties of standard symmetries for what concerns 
reduction of a scalar ODE [12]. 
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In the present note, we extend the concept of A-symmetries to the case of 
partial differential equations (PDEs), and to systems. In order to obtain such 
an extension, we found it convenient to characterize A-prolongations in J^"'M, 
where (M, tt, B) is the space of dependent and independent variables seen as a 
bundle over the space B of independent variables, in a geometrical way; once 
this characterization is obtained, it is promptly extended from the ODEs case 
of i? = R to the PDEs case B = IV, and to systems. 

In the scalar PDE case the transformations of interest depend on a semibasic 
one form fi ~ A^dx*, the functions Aj being such to satisfy a compatibility con- 
dition. The transformations of this class leaving invariant the solution manifold 
for an equation A will be said to be /x-symmetries, or A^-symmetries, of A. In 
the case of systems, the form /x takes value in a Lie algebra. 

We will thus be able to obtain a sound definition of /x-prolongations and 
/x-symmetries of (systems of) PDEs. We will also show that, in analogy with 
the ODE case, /x-symmetries are as useful as standard symmetries in what 
concerns the symmetry reduction, and the determination of invariant solutions, 
of (systems of) PDEs. Our approach will suffer from the same limitations as 
the standard PDE symmetry reduction method. 
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1 Prolongations and contact structure 

In this section we fix our notation (mainly following the usual one in the field, 
see e.g. [10]) and recall some basic definition [4, 10, 13, 15]. 

Let us consider a space M = B x U with coordinates x G B ~ IW and 
u € U c^ R*?; when setting a differential equation in this space, we will think of 
the X as independent variables, and the u as dependent ones. 

Thus, more precisely, M will be the total space of a (trivial) linear bundle 
(M, TT, B) over the base space i?, with fiber tt~^{x) = U . 

Given a bundle P, we will denote T[P] the set of sections of this bundle, and 
by X[P] the set of vector fields in P. 

The bundle M can be pro^onf^ed to the fc-th jet bundle {j'^^^ M,iTk,B), with 
j'^^^M = M; the total space of the jet bundle is also called the jet space, for 
short. 

The jet space J^'^^M is naturally equipped with a canonical contact struc- 
ture £, i.e. the module generated by the set of canonical contact one-forms 

rj := du} - ul^dx- 

with a = l,...,g, \J\ = 0, ...,k — 1. 



The contact structure in J^^^ M defines a field of (p + g)-dimensional linear 
spaces in J'^^^M C T((J'^'^~^^M), the contact distribution, corresponding to the 
tangent subspace spanned by vector fields Y £ A'[j('^^M] annihilated by the 
contact forms, i.e. such that Y _\9 = for any contact form 9. The general 
form of such vector fields is, as well known, Y = Y^ C^l + ^• 

Here Di is the total derivative [4, 10, 13, 15] with respect to a;*, D^ ' its 
truncation to the fc-th jet space, and t^ is a generic vector field in X[j'-'^^M], 
vertical for the fibration nk.k-i '■ J^'^^M -^ J^'^^^^M (the latter will not appear 
if we work with infinite-order prolongations; it will however disappear when we 
deal with a given differential equations and symmetry vector fields for it). The 

(k) 

operator D^ reads 

q fe-1 
a=l |J|=1 

In the following we will write, for ease of notation, simply Di instead of Dl . 
A vector field X e X[M] can be written, in the {x, u) coordinates, as 

This can be uniquely prolonged to a vector field X^ '^ in J^^'M by requiring it 
preserves the contact structure (the precise meaning of this will be defined in 
a moment). The prolongation formula [4, 10, 13, 15] is indeed expressing this 
condition. 

We write a vector field in j'-'^^M as 

F = X + V ^".-^ 

where X is as above, J = ji,...,jp is a multiindex, and the sum is over all 
multiindices of modulus \J\ = ji + ... + jp up to the order of the jet space. We 
also write Dj for the total derivative Z)"'\...Z)^p, and Uj for Dju'^; moreover 
uj^i will denote DiUj. 

Then Y is the prolongation of X if and only if the coefficients ^ } satisfy the 
prolongation formula 

*} = Dj^'^ - Dj{CO + eOjuf . (1) 

This is also recast in recursive form. We denote by J = J + e^ the multiindex 
with entries ji = ji+6ik, and for short uj^k := uj+e^, *},fc := ^J+e^- Then (1) 
is equivalent to 

*lfe = Dk^j - ul^Dkr (2) 

with *g = ^" (see e.g. sect. 2.3 of [10]). 



Let us now discuss the geometrical aspects of the prolongation operation 
in terms of contact structures [4, 7, 10, 11, 13]; these will be useful for our 
subsequent generalization. 

Note prcliminarly that for any function / : J^^^^> M ^ R we can write 

d/ = (A/)dx* + d[f] (3) 

where i?[/] G £ is some contact form whose explicit expression (easy to compute) 
is irrelevant here. 

Definition 1. Let F be a vector field on J^'^^M. We say that Y preserves 
the contact structure if Cy '■ £ ^ £■ ♦ 

Proposition 1. The vector field Y G X[J^''^M], projecting to a vector field 
X G X[M] on M , is the prolongation of a vector field X G X[M] if and only if 
it preserves the contact structure in J^'^' M . 

Proof. This is a classical result, see e.g. [7, 10, 11, 13]. <() 

Lemma 1. The vector field Y preserves the contact structure £ if and only if, 
for any -d ^ £ and any i = I, ...p, ([-Di, Y\) _!'(? = 0. 

Proof. Write d°j and Y as above, and note that Di ~ di + Uj,-{d/duj), where 
of course di = d/dx^. With this notation and standard computations, 

[A,r] = {D^Ddm + (A*"; - ^%){d/du'}) ■ (4) 

hence we get [A,i^]-I^} = -^},, + (A*} - w},™(AC")), which vanishes if 
and only if the 4"" satisfy the recursive prolongation formula (2). () 

Corollary. The vector field Y preserves the contact structure £ if and only if 
[Di,Y\ = K^Dm + V for some /i™ G A°{J^'''>M) and V a vertical vector field 
for the fibration iTkM-i ■ J^'^'^M -^ J^''~^'>M. 

Proof. The vector fields Dm span the set of non- vertical vector fields (for the 
fibration T:k,k-i) in the annihilator of the contact forms. Alternatively, this 
follows at once from (3), with h™ = Di^™. <C> 

2 ODEs: A-prolongations and A-symmetries 

In this section we will restrict to the case of scalar ODEs, i.e. to the case where 
the bundle {M, n, B) has _B = R as base space and 7r~^(&) = R as fiber. We will 
characterize in geometrical terms, i.e. in terms of their action on the contact 
structure, the A-prolongations introduced by Muriel and Romero [8] (see also 
[9] and [5]), and further studied by Pucci and Saccomandi [12]. 

We simply write u„ for -D"u, and similarly for ^„. The standard contact 
forms in J'-'^'M will be i?„ = du„ — Un+idx, with n = 0, ..., fc — 1. 

We start by recalling the definition of A-prolongations and A-symmetries as 
given by Muriel and Romero, using an obvious notation for x-derivatives of the 
u. 



Definition 2. Let X ~ £_{d/dx) + (p(d/du) be a vector field on M, and Y = 
X + J2n=i *«(9/(9u„) a vector field on j'-'^^M. Let A : j'-^^M ^ R be a smooth 
function. We say that Y is the \- prolongation of X if its coefficients satisfy the 
A-prolongation formula 

*„+i = [(i?. + A)*„] -M„+i [(/?, + AK] (5) 

for all n — Q, ..., fc — 1. Jit 

Definition 3. Let A be a fc-th order ODE for u = u{x), u E U = H, and 
let (M = U X B,Tr,B) be the corresponding variables bundle. Let the vector 
field Y in J^'^^M be the A-prolongation of the Lie-point vector field X in M. 
Then we say that X is a A-symmetry of A if and only if Y is tangent to the 
solution manifold ^a, i.e. iff there is a smooth function $ on J^^ M such that 

r(A) = $A. 4k 

Remark 1. We stress that in this note we take A : J^^^M -^ R, which guaran- 
tees that the A prolongation of a Lie-point vector field in M is a proper vector 
field in each j'^'^^M. One could also consider A : j'^'^^M -^ R, obtaining obvious 
generalizations of the results given here. In this case the A-prolongations of X 
would be generalized vector fields in each j'^'^^M with n > even if X is a 
Lie-point vector field. The same applies to the /i-prolongations to be considered 
in later sections. 

We will not discuss here the relevance of A-symmetries, referring to [8, 12]; 
we just recall that they are as useful as standard ones in that one can perform 
symmetry reduction to the same extent as for standard symmetries [8] . 

The basic property of A-prolongations behind this feature was clearly pointed 
out by Pucci and Saccomandi [12], and can be expressed in terms of the char- 
acteristics of the vector fields Y which are A-prolongations of X . 

Given the vector bundle {M,tt,B), we choose a distinguished smooth real 
function X{x,u,Ux) '■ J'^^'M -^ R. We note for later discussion that to this 
is associated a semibasic one-form fx e A^(J(^)M), i.e. the one-form /i = 
A(x, u, Ux)dx 

Definition 4. Let F be a vector field on the contact manifold (J'^'^^M, 5), and 
A e A"(J(^)M) a smooth function on M. We say that Y A-preserves the 
contact structure if, for any contact one- form 9 G £, 

£y{0) + (y_i6')Adx = e (6) 

for some contact one- form 9 E £. 4^ 

Theorem 1. Let (M, tt, _B) be a bundle over the real line B = R witti fiber 
7r~^(x) = R, and let £ be the standard contact structure in J^^'M . Let Y be 
a vector field on the jet space J^^'M , which projects to a vector field X on 



M . Then Y is the X-prolongation of X if and only if it X-preserves the contact 
structure. 

Proof. We write a general vector field on j'^^^M asY ^ C^a;+X]m=o ^ ^.{9 / dum)] 
as the contact forms are ■!?„ = dM„ — M„+idx {n = 0, ..., fc— 1), we have by explicit 
computation 

with 9 a contact form. Thus (6) is satisfied if and only if the ^„ satisfy the 
A-prolongation formula (5). <C> 

We can also provide an alternative characterization of A-prolonged vector 
fields, similarly to what we did for standard prolongations in lemma 1. 

Lemma 2. Let Y he a vector field on the jet space J^'^^M, with {M,tt,B) 
a vector bundle over the real line B = R, and let £ he the standard contact 
structure in J^^' M . Then Y is the \-prolongation of a vector field X on M if 
and only if, for any d £ E, 

[D,,Y]_^d = \{Y_^d) . (7) 



Proof. Looking at the proof of lemma 1, [D^, Y]_\d is given by (4) specialized 
to the case p = \: with the obvious notation m„ :— D"u (and similarly for 
\l/„) we have [£'a;,y] = — ^„+i + (Dx'i'n — Un+iD^S,); on the other hand, it 
is easy to check that y_l??„ = 'i'n — Un+iS,- Thus cq. (7) is equivalent to 
\l/„+i = [{Dx + A)^„] — Un+i[{Dx + X)S,], i.e. to the A-prolongation formula (5). 



Corollary. In the hypotheses of lemma 2, Y is the X-prolongation of a vector 
field X on M , if and only if[Dx, Y] = XY+hD^+V with A, h scalar functions on 
J^^^M and V a vertical vector field for the fibration TTk,k-i ■ J^'^'^M -^ J^'^^^^ M . 

Remark 2. Theorem 1 shows that our geometrical formulation, i.e. definition 4, 
is equivalent to the standard (analytical) one, i.e. definition 2. The advantage 
of our formulation is twofold: we have a geometrical characterization of A- 
prolongations (A-symmetries), and moreover this is readily extended from the 
ODEs to the PDEs case. As we discuss later on, we can moreover generalize the 
standard symmetry reduction method for PDEs to an analogous A-symmetry 
reduction. We will also show that this definition and the reduction procedure 
extend to systems of PDEs. © 

3 PDEs: yU-prolongations and ^u- symmetries 

In this section we extend our approach to A-prolongations and A-symmetries to 
the case of scalar PDEs {p independent variables); the case of PDE systems will 
be dealt with in section 5 below. 



The role of the scalar function A will now be played by an array of p smooth 
functions \i : j'^^^M — > R (remark 1 holds also in this context), which will be 
the components of a semibasic form ^ S h}{J^^^M). The only additional ingre- 
dient required in the multi-dimensional (PDE) case is a compatibility condition 
between the semibasic form /x and the contact structure - this is eq.(lO) below 
- automatically satisfied in the ODE case. 

Actually, our formulation of A-prolongations in the ODE case was such that 
the results, and even their proofs, are the same also in the PDE case - except 
of course for the appearance of new variables. 

In view of our geometric approach it is convenient to focus on the form ^ 
rather than on the g-ple of smooth functions Xi . We will thus call the analogue 
of A-prolongations and A-symmetries in the PDE frame, /x-prolongations and 
/K-symmetries. 

We equip ( j'^-'M, tt, B) with a distinguished semibasic one- form /i, 

H = \,dx' . (8) 

We require that /u is compatible with the contact structure defined in J'^'^^M, for 
fc > 2, in the sense that 

All G J{£) , (9) 

where J{S) is the Cartan ideal generated by £ (we recall that a two-form a is 
in J{£) if and only if a = p'^ A 'dj for some one-forms p^). 

It should be noted that this condition does not appear when we deal with 
first order equations, i.e. with first order /i-prolongations. We note also that 
for p = 1 eq.(9) is automatically satisfied: indeed, d/i = {d\/du)du A dx + 
{dX/du,j.)dux f\dx = {d\/du)-dQ A dx -I- {d\/du,j.)i)i A dx. 

Lemma 3. Condition (9) is equivalent to 

DiXj - DjXi = . (10) 

This is in turn equivalent to the condition that the operators V^ := Di + Xi 
commute, [Vi, Vj] = 0. 

Proof. As Xi is a function on j'^'^'' M , we have dp — {dXj / dx^)dx^ A dx-' -|- 
{dXj / du)du A dx^ + (dXj/dui)dui A dx-', i.e. 

d^ — [{dXjIdx^) + Ui{dXj/du) + Uik{dXj/duk)]dx^ A dx-' + 
+ {dXj/du)i9o A dx^' + {dXj/dui)d^ A dx^' . 

The two latter terms are of course in J{£), while no form dx* A dx-' belongs 
to J{£)\ thus, (9) is satisfied if and only if the coefficients of all these terms 
vanish. This condition is precisely (10). (Note this extends to the case where 
p is semibasic for (j("\7r„,i3), see remark 1.) The equivalence of this with 
[Vi, Vj] = follows from the definition of V^. () 

Definition 5. Let F be a vector field on the contact manifold {J^'^' M,£), and 
p a semibasic form on M compatible with £. We say that Y /i-preserves the 



contact structure if, for any 6 E £, there is a form E £ such that 

CriO) + iY^e)fi = 9. (11) 

Definition 6. A vector field Y in J^'^^M which projects to X in M and which 
/x-preserves the contact structure is said to be the /i-prolongation of order k, 
or the fc-th /x-prolongation, of X. Jit 

Theorem 2. Let Y be a vector field on the jet space j'^^^ M , with {M,t:,B) a 
vector bundle over B = R^, written in coordinates as 

with X = S^^(d/dx^) + ip{d/du) a vector field on M . Let £ be the standard contact 
structure in J^ ' M , and ji = A^da;* a semibasic one-form on (J^^' MjTT, B), com- 
patible with £ . Then Y is the ji-prolongation of X if and only if its coefficients 
(with ^0 = v) satisfy the /i-prolongation formula 

^j,, - (A + A,)*j - wj,„ (A + K)r ■ (12) 



Proof. The standard contact forms in J^^^M are -Dj = duj — uj^dx^, with 
\J\ = 0, ..., fc — 1. Thus, as already computed in the proof of proposition (1), 
'Cy(^j) == (-*j,i + A*j - u,7,„iAC")da;* + 6 with 6 a contact form. On the 
other hand, it is easy to compute that y _l??j = 'i'j — uj.mC™- Therefore, 

This is a contact form if and only if the coefficients of all the da;* vanish, i.e. if 
and only if (12) is satisfied. ^ 

Remark 3. Condition (9) arises from the following: consider the multiindices 
J — (ji, ■■■,jp) and L = {£i, ...J.p) with tg = js + Si^s + Sk,s', the coefficient ^^ 
can be obtained from '^j by applying twice formula (12), but we can proceed 
in two different ways, i.e. pass first from ^Pj to '^ja and then to '^l, or pass 
first from ^j to '^j,k and then to '^l. Needless to say, the result must be the 
same in the two cases, and this is the compatibility condition for the A^. By 
explicit computation this is just (10), equivalent to (9) by lemma 3. 

As for standard and A-prolongations, /i-prolongations have a specific be- 
haviour for what concerns their commutation with the total derivatives Di. 

Lemma 4. // Y is the pL-prolongation of a Lie-point vector field X , with /i = 
Xidx^ , then for any contact form d, 

[D„Y]_^d = A,(yj7?) . (13) 



Proof. In the proof of lemma 1 we have computed [Di,Y]_\i!)j = — ^j^i + 
Di'i/j — uj^mDiS."^', needless to say, Y _\i}j = *j — uj^mS,"^ and thus (13) is 
equivalent to the /^-prolongation formula (12). (} 

Corollary. In the hypotheses of lemma 4, Y is the ^-prolongation of a vector 
field X on M , if and only if [D^, Y] ~ X^Y + h^Dm + V with Xi, K^ scalar 
functions on J^^' M and V a vertical vector field for the fibration of J^^'M over 

It is quite remarkable that a simple relation exists between the /i-prolongation 
of a vector field and its ordinary prolongation. In order to discuss this relation, 
we write X = ^"^{d/dx^) + ip{d/du) for the vector field in M, and denote its 
ordinary prolongations as X^'''' = X ^"^ j{d/duj), while its /i-prolongations are 
denoted a.sY = X + '^ j{d/duj). The form /i is written, as usual, /i = A^da;*, 
and of course '^ j ^ "^ j when all the Xi (or at least all those for i such that 
ji ^ 0) vanish. 

As well known [4, 10, 13, 15], the vector field X can be cast in evolutionary 
form as Xq := Q{d/du), with Q := ip ~ Ui^"^ . 

The equations DjQ — 0, with \J\ = 0,...,fc — 1 identify the X-invariant 
space Ix C J'^^^M. We denote by T the module over C°°{J^'^^M) generated by 
the DjQ, i.e. the set of functions F which can be written as F = c^ DjQ for 
some smooth functions c^ : J^^^ M -^ R, and by ^(™) <Z T = JFC^) those which 
depend only on variables (x,!*'™)), m < k. Needless to say, Di : jf(™~i) -^ 

Theorem 3. Let X,Y,^ be as above. Write *j as *j = vj/j + Fj. Then the 
functions Fj satisfy the recursion relation (with Fq = 0) 

FJ,^ = (A + K)Fj + KDjQ . (14) 



Proof. In order to show that the statement of the theorem holds at all orders, 
we proceed recursively: we suppose (14) holds for all \J\ < h, and wish to prove 
that it holds also for \J\ = h. Any J of order h can be written as J + e^ for 
some i and some J of order h — 1; formula (14) holds for '^j. Thus, by the 
/^-prolongation formula, 

*j,^ = (A + \)^.i - wxm(A + \)e' + (A + A,)Fj 

= [A*,/ - W J,™ AC"] + A^ [*,/ - Uj^U"'] + {D^ + Xi)Fj . 

The first term is just '^ j^i, and the last is already in the form appearing in (14); 
so we have to look only at the second one. 

Take an s such that js 7^ in the multiindex J, and write K = J — e^ . Then, 
using the standard prolongation formula. 



We can then repeat the procedure on any index q such that K = J — Cg has a 
nonzero q entry, and so on. In the end, recalhng that \&o = ^, we have 

*,/ - Uj,,nr = Dj (^ - UrnC) = DjQ ', 

this also follows from the formula for prolongation of the evolutionary represen- 
tative of X. Going back to our computation, we have thus shown that 

*J,» = ^.U + \^DJQ + (A + \^)FJ . 

This shows that if (14) is satisfied at order ft, — 1, it is also satisfied at order h. 
It is easy to check that (14) holds at order one, i.e. for \J\ = 0: indeed, by 
the /Lt-prolongation formula (12) and the ordinary prolongation formula (1), 

*, = {D, + \,)^-u^{D, + \i)^"' 

= (A(^ - u™Ae") + A, (</? - u„^") 

We conclude that (14) holds at all orders. () 

This theorem provides an economic way of computing /i-prolongations of X 
if we already know its ordinary prolongations. Theorem 3 also has a rather 
obvious consequence, which will be relevant in the following. 

Lemma 5. Let X be a vector field on M , £ the standard contact structure 
on J^^^'M , and /i any semibasic form on M compatible with £. Then: (i) 
the ^-prolongation Y of X coincides with the ordinary prolongation X^''' on 
the invariant space Tx; (ii) the space Tx C J^^'M is invariant under the fi- 
prolongations of X , for any semibasic form /i compatible with £. 

Proof. By definitions, any function F ^ T vanishes identically on Ix- Thus 
(14) guarantees that \I'j = ^j on Xx^ i.e. proves point (i). 

As for (ii)^ this is a known property of standard prolongations, easily checked 
by using the evolutionary representative of X, Xq :— Q{d/du). Its prolongation 

is Xq = {DjQ){d/duj), where the sum is over all multiindices with | J| < A:, 

and XC^) = X^Q^ + f A- Thus, X^^) reduces to I^ = f Di on Ix] and W is 
obviously tangent to Ix ■ ^ 

Finally, we define /i-symmetrics of a PDE as Lie-point vector fields whose 
^-prolongation is a symmetry of the equation. 

Definition 7. Let X be a vector field on M, and let Y e A'[J('=)M] be its 
/i-prolongation of order k. Let A be a differential equation of order k in M, 
A ■- F(a;, «(*=)) = 0, and S C J(*=)M be the solution manifold for A. If 
Y : S ^ TS, we say that X is a /i-symmetry for A. If Y leaves invariant each 
level manifold for F, we say that X is a strong /i-symmetry for A. ^ 

Remark 4. Note that if we look for /i-symmetries of a given equation A, we 
can accept forms /i which do not satisfy (9) on the whole jet space J^'^'M, but 
only on the solution submanifold ^a C J^^-^M. 
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Remark 5. Given a form fi = A^dx*, we consider exponential vector fields 

X ^ e/ ^ • Xo 

where Xq is a vector field on M; note that ii fi — {DiP)dx^, in which case (10) 
is automatically satisfied, then X — c^Xq; in general X is a formal expression. 
For a general fi, consider an equation A such that (10) is satisfied on Sa, see 
the remark above. Then we have the following result (see [2] for a proof and 
extensions): X is a (in general, nonlocal) symmetry for A if and only if Xq is a 
/x-symmctry for A. This extends a result by Muriel and Romero [8]. 

The relevant point is that /i-symmetries can be used to obtain group- 
invariant solutions, i.e. one can introduce /i-symmetry reductions of PDEs 
and obtain invariant solutions to the original PDE from these, by the same 
method as for standard symmetries. 

Note that in this way we parallel again the ODE case, where it was proven 
by Muriel and Romero and by Pucci and Saccomandi that A-symmetries are as 
good as standard ones for reduction of the equation. 

4 The fi-syuiuietry reduction method for PDEs 

As well known, symmetry reduction for PDEs is conceptually different from 
symmetry reduction for ODEs: while in the latter case it yields a reduced 
equation whose solutions provide, together with an integration, the most general 
solution to the original ODE, in the PDE case the reduced equation provides 
only the symmetry-invariant solutions to the original PDE. 

4.1 The PDE reduction method 

In this subsection we briefly recall (using the notation introduced so far) sym- 
metry reduction for scalar PDEs in the case of standard symmetries; this is 
discussed in detail in a number of textbooks and research papers, see e.g. 
[4, 10, 13, 15]. We will just discuss reduction under a single vector field, rather 
than a general (i.e. higher dimensional) Lie algebra. 

Consider a PDE of order k A, which we may think in the form F{x, u*-*^' ) — 
with F : j'^^^ M ^ R a smooth scalar function. Let the Lie-point vector field 
X in M, with prolongation X'^^'^ in J^'^^M, be a (standard) symmetry for A. 
Then we proceed as follows, following Olver. (For more details, see e.g. the 
discussion in chapter 3 of [10]). 

First of all we pass to symmetry-adapted coordinates in M. In practice, 
we have to determine a set of p independent invariants for X in M, which we 
will denote as {y^, ...,yP^^,v): these will be our X-invariant coordinates, and 
essentially identify the G-orbits, while the remaining coordinate a will be acted 
upon by G. In other words, G-orbits will correspond to fixed value of (y, v) 
coordinates and to a taking values in a certain subset of the real line (thus 
{y,v) are coordinates on the orbit space il — M/G, see [10]) 
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The invariants will be given by some functions y* — "rfix^ u) {i = 1, ...,p — 1) 
and V — C(x, u) of x^, ..., x^ and u. If X acts transversally, we can invert these 
for X and u as functions of (y, w;cr), i.e. write a;* = x*(y,w;(T) {i ~ 1, ...,p) and 
u = I3{y,v;a). 

If now we decide to see the (y; a) as independent variables and the v as the 
dependent one, we can use the chain rule to express x-derivatives of m as ct and 
y-derivatives of v"^ . Using these, we can finally write A in terms of the (y, u; a) 
coordinates and derivatives of v in the y and a; this will turn out to be an 
equation which, when subject to the side condition dv/da = 0, is independent 
of a. The condition dv/da = expresses the fact that the solutions are required 
to be invariant under X, i.e. the equation obtained in this way represents 
the restriction of A to the space of G-invariant functions, and therefore it is 
sometimes also denoted as A/G. 

Suppose we are able to determine some solution v — ^(y) to the reduced 
equation; we can write this in terms of the {x,u) coordinates as (^{x,u) = 
^[r]{x,u)], which yields implicitly u = f{x): this is the corresponding X- 
invariant solution to the original equation A in the original coordinates. 

Remark 6. The symmetry reduction method for PDEs can also be seen in a 
slightly different way: if we look for X-invariant solution u = ,f{x) to A, we 
determine the characteristic Q — (p — Ui£/^ of the vector field X, and supplement 
A with the equations Ej := DjQ = with \J\ = 0, ..., fc — 1. The equation E^ 
requires that the evolutionary representative Xq — Q{d/du) vanish on 7/, i.e. 
that u is X-invariant, and all the equations with | J| > are just differential 
consequences of this. The X-invariant solutions to A are in one to one corre- 
spondence with the solutions to the system A(x) := {A;_Ej}. See e.g. [15] for 
details, and for how this approach is used in a more general context. 

Remark 7. We stress that the standard method discussed here applies un- 
der a nondegeneracy (transversality) condition, guaranteeing a certain Jacobian 
admits an inverse. When this is not the case - as it happens in a number of 
physically relevant cases - the treatment should go through the approach de- 
veloped by Anderson, Fels and Torre [1]. See also [6] for the case of partial 
transversality. We also stress that this method is justified only if the (possibly, 
only local) one-parameter group G generated by X has regular action in M, i.e. 
the G-orbits are regular embedded submanifolds of M [3] . In the following we 
tacitly assume both conditions mentioned here are satisfied. 

4.2 On the justification of the method 

The method described above is rigorously justified in chapter 3 of [10], to which 
we refer for details. In this subsection we just recall what is the key step in the 

^It is maybe worth recalling that this computation can be described also in a slightly 
different, but equivalent, way: that is, we write du = Uidx^ on the one hand, and du = 
d[f3{y,v;a)\ on the other. We then expand the latter as du = f3jdy^ + Pvdv + fS^dcr, where 
of course /3j = d/S/dy^ , and substitute for dv as dv = Vjdy^ + v„-dcr. Comparing the two 
expressions for du, we obtain the expression for Ui in terms of Vj and Do-- 
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proof, as we will have to prove a similar property also holds for /i-prolongations 
in order to justify the extension of this method to /^-symmetries. 

We recall that a function u = f{x^,...,xf) corresponds to a section 7/ G 
r[M], the space of sections for the bundle {M,tt,B), i.e. jf — {{x,u) : u = 
fix)}. This is uniquely prolonged to a section ^V' G r[j'^'^)M]; ^^ ' is the 
unique lift of the curve 7/ in M to a curve in J^^> Ad which (i) projects down to 
7/ in M, and [ii) is everywhere tangent to the field of contact linear spaces. 

When we act with a vector field X = £^'^{d/dx^) + (p{d/du) on M, at the 
infinitesimal level, the section 7/ is mapped into 7-7 with 

/(x) == f{x) + e[tp{x,u)-£_\x,u)dif{x)]u=f(a:) + o{e) . 

Thus a function u = f{x) is invariant under the action of the vector field X 
in M if and only if Q{x) := Q[x, f{x)] = (with Q the characteristic of X). 

We consider the equation Eq := Q — Q and all of its differential consequences 
Ej :— DjQ — for \J\ < k; this identifies the invariant manifold Xx C j'^^^M. 
Passing to the evolutionary representative Xq = Q{d/du) of X, it is obvious 
that Xq and its prolongations vanish on Xx- The X-invariant solutions to 
A will be the solutions to the system A(x) made of A and of the invariance 
condition: 

F(x,wW) = , . 

DjQ = (|J|=0,...,fc-l) . ^ ^ 

We denote the solution manifold to this system as Sx C 2x C J^'^^M. The 
invariance ofSx, as discussed by Olver [10], guarantees that the method recalled 
above is justified. 

Recall now that the prolongations of X and Xq satisfy 

X^'^) = Xf ^ + f i?f ^ . (16) 

Lemma 6. The (standard) prolongation X^'^' of X reduces to £,^Di onXx, and 
is tangent to Sx- 

(k) 
Proof. The field Xq ' vanishes on 2x because of the equations Ej, and the Dj 

are symmetries of any system, as the differential consequences of any equation 

of the system are satisfied by solutions to the system. By (16), this proves the 

claim. ^ 

4.3 Reduction of PDEs under yu-symmetries 

In the case of /i-symmetries of PDEs, we can proceed exactly in the same way 
as for standard symmetries in order to determine G-invariant solutions. 

Note that the step consisting in the introduction of symmetry-adapted co- 
ordinates is exactly the same; the difference lies of course in the step connected 
to the prolongation structure. 
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We describe here how the standard symmetry reduction method is formu- 
lated to deal with ^-symmetries. We suppose that X is a /x-symmetry of A, 
acting transversally for the fibration (M, tt, B), and denote the /^-prolongation 
of X as ye X[J'^^^M]. 

First of all we pass to symmetry-adapted coordinates (y, u; a) in M, as in the 
standard case. We retain the notation introduced in subsection 1. We further 
proceed as there, i.e. use the chain rule to express x-derivatives of the m as cr 
and y-derivatives of the v. Using these, we can finally write A in terms of the 
{y,v;a) coordinates and their derivatives. 

Again, looking for X-invariant solutions means supplementing the equation 
with the side condition dv/da = 0, or with the conditions DjQ = 0, see eq.(15) 
above, in the original coordinates. 

Now the point is that if the equation thus obtained is independent of a, we 
have indeed obtained a symmetry reduction of the original equation. In this 
case solutions v ~ $(y) to the reduced equation can be written in terms of 
the {x,u) coordinates as (^{x,u) — ^[ri{x,u)] and yield implicitly u — /(x), the 
corresponding X-invariant solution to the original equation. 

However, the vector field Y is not the ordinary prolongation of X , and thus 
we are not apriori guaranteed it leaves Sx or Tx invariant. Thus, in order to 
justify the method sketched above - i.e. in order to prove that the standard 
PDE reduction method still applies in the case of /i-symmetries - we have to 
prove the following theorem 4. Note that the only difference with respect to 
the standard case will be that it is the vector field Y, and not the ordinary 
prolongation X^'^^ oi X , to be tangent to the solution manifold of A in j'-'^^M. 

Theorem 4. Let A be a scalar PDE of order k for u = u{x^ , ..., x^). Let X = 
£,''{d/dx'^) + (p{d/du) be a vector field on M, with characteristic Q := ip — Ui^"^ , 
and let Y be the fi-prolongation of order k of X . If X is a ji- symmetry for A, 
then Y : Sx -^ TSx, where Sx C J^^'M is the solution manifold for the system 
/S.X made of A and of Ej :— DjQ — for all J with \J\ = 0, ...,k — 1. 

Proof. Recall that Sx is the intersection of the solution manifold Sq to A with 
the X-invariant set Xx (see remark 3 above, or [15]). The former is y-invariant 
by assumption, as X is a /x-symmetry of A; the y-invariance oiXx is guaranteed 
by lemma 5 above. Therefore the proof for the standard case [10] extends to 
the present setting. (} 

Remark 8. The property Y : Tx -^ TXx can be shown in a alternative 
way without resorting to comparison with the standard case, i.e. using the 
geometrical characterization of /x-prolonged vector fields, as follows. 

Denote byXj^^ C J^'^^ M the set of points identified by Ej for |J| < m. We 
first show that iiX^ is invariant under Y, then Xj^ is also K-invariant (for 
m — 0, ...,k — 2). Note that F-invariance oi 2^ means that for all |J| < m 
there are functions P^ such that Y{DjQ) — J2\k\=o/^^ ^^Q- 

We have Y[D^{D.jQ)] = [Y, Di]{DjQ) - D^{Y{DjQ)); from the corollary to 
lemma 4 this reads KY{DjQ) + hlDs{DjQ) - D^{Y{DjQ)) + V(DjQ), with 
V — J2\K\=k^^i^/'^^K)- The first term is in Xj^"-^ by hypothesis, while the 
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second and third ones are by definition in X)^ ' . The last term vanishes since 
DjQ does not contain u derivatives of order greater than m + 1, and m < k — 2. 

The proof of K-invariance of Xx is hence reduced to proving Y-invariance 
of Xjj- , i.e. of the manifold identified by Q = 0; as for X a Lie-point vector 
field Q depends only on first order derivatives, it suffices to consider the first 
/i-prolongation of X, which is just X^^^ + XiQdu^- It is well known that Q = 
is invariant under the ordinary prolongation X^^' , and of course the other term 
vanishes on Q = 0. 

This proves y-invariancc of Xj^- and hence, by the recursive argument given 
above, of all the Xj^ with m = 0,1, ...,k — 1. 

The recursive property considered here can be seen as a counterpart in the 
PDE case to the recursive property discussed by Pucci and Saccomandi as char- 
acterizing the A-prolongations as telescopic vector fields in the ODE case [2]. 




5 Systems of PDEs 

In this section we extend ^-prolongations to the case of g > 1 dependent vari- 
ables. We will assume that the dependent variables u take value in the vector 
space U = R«, and B = Rp. 

It will be natural in this context to consider differential forms taking values 
in the space G = 5^((?), the Lie algebra of the group G = GL{q). Thus we 
will deal with matrix-valued differential forms (or more generally Lie-algebra 
valued differential forms), see e.g. [14]. The form jjl will now be written in local 
coordinates as 

^l := {K^)lAx' (17) 

where A^ : J^^' M -^ Q arc smooth g-dimcnsional real matrix functions. Note 
that remark 1 applies also to this case. 

5.1 yU-prolongations in vector framework 

In this case we generalize condition (6) to the following (19); we will then define 
^-prolongation in the same way as in the scalar case. 

In the vector case, we see the contact structure Q (we use a different sym- 
bol than in the scalar case to emphasize we deal with vector-valued forms) as 
spanned by vector-valued one-forms d,j = (i?j, ..., z?j) e R'' (E) A-'^(M), where 

rj = du} - ul^dx"^ . (18) 

Definition 5'. We say that Y ^-preserves tlie contact structure Q, with 
^ given by (17), if for any vector- valued contact forms ^ G O, there is a vector- 
valued contact forms i3 e Q such that 

/:y(r) +(y_i [(A.)fc/]) dx' = r . (19) 
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4k 

Definition 6'. A vector field Y in J^^^M wliich projects to X in M and which 
^-preserves the contact structure is said to be the /i-prolongation of order k, 
or the fc-th ^-prolongation, of X. Jit 

In order to discuss vector fields in J^^' M which are /i-prolongations of vector 
fields in M, it will be convenient to agree on a general notation. That is, we 
write a general vector field in j'^^^M in the form 

^ - ^^^ + *"^4 • ^'°^ 

Theorem 5. The vector field Y fi-preserves the standard contact structure <d 
if and only if its coefficients satisfy the vector /i-prolongation formula 

*},. = KD, + (A,)?] ^'j - u\^ [5tD, + {K)l] e . (21) 



Proof. This follows easily by a computation analogous to that in the proof of 
theorem 1. <^ 

Theorem 6. Let X = £^^{d/dx'^) + ip"-(d/du"-) be a vector field in M. Let 
jjL = (Ai)^dx* he a Q-valued semibasic one-form. Then the coefficients ^"j of 
the ^-prolongation Y of X are expressed in terms of the coefficients ^ j of the 
ordinary prolongation of the same vector field X as ^ " = 4'" + _Fj where the 
difference terms Fj satisfy the recursion relation (with F(f = Q) 

Fh = KA + (A.)?] F^j + {k,)tDjQ' . (22) 



Proof. Follow the scheme used in the proof of theorem 3. ■O 

Similarly to what happens for the A^ in the scalar case, see remark 3, the 
(matrix) coefficients A^ of the form p, are not completely arbitrary, as they must 
satisfy some compatibility condition. It is convenient to define the (matrix) 
operators V^ := IDi + A^. 

Theorem 7. The compatibility condition for the matrix coefficients A^ of the 
Q-valued form /i — (Ai)^dx* reads 

A A, - DjK, + [A,, A,] ^ (23) 

for all i,j — 1, ...,p. This is equivalent to [ V^ , Vj ] ~ 0. 

Proof. We will use the shorthand notation (Vj)JJ = [S'j^D, + (Aj)JJ]. With this, 
the vector /^-prolongation formula (21) reads 

*},fc - (^k)ttj - 4™(Vfc)?r ; 
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applying this twice, we get 

^Im - [(v.v,)gvi/';_4^(v,v,)^r] - K.,™(v,)g + 4fe,„(v,)^]r ; 

note the second square bracket is symmetric in the indices i, k. Thus 

As for the commutator [Vi, Vfe], this is easily computed to be 
[V„Vfc] == AAfc - DkA, + [A,,Afe] , 

i.e. the expression given in the statement, see (23). {)> 

Remark 9. One could consider matrices A^ belonging to a gauged Lie alge- 
bra. By this we mean that A^ ~ Af (x,m^^)) ife, with Ai : J^^^^M -^ R smooth 
functions and where the Lk {k — 1, ...,r) are generators of a (matrix) Lie alge- 
bra G, so that [Li,Lj] = c'^jLk- In this case the compatibility condition reads 
[{DiX'j — DjX^) + c^fcAfA^] Lk = 0; the term in square brackets must vanish 
for each k. 

5.2 yU-symmetries and reduction of PDE systems 

We define /u-symmetries as in the scalar case; that is, A is a ^-symmetry of a 
given PDEs system if its /i-prolongation is tangent to the solution manifold of 
the system. For scalar equations, /i-symmetries can be used to obtain invariant 
solutions; the same holds for the vector case. 

We will assume without further mention that X satisfy the transversality 
condition in the bundle (M, tt, i?) [10]; see remark 7. 

Definition 7'. Let {M,tt,B) be a vector bundle over the p-dimensional man- 
ifold B, with fiber 7r-i(x) = ^7 = R"*. Let A = {Ai,...,Ar} be a system of 
PDEs of order n for u" ~ m"(x), a — 1, ...,9, x — (x^, ...,2;'') G B, with solution 
manifold Sa C J^^'M. Let A be a vector field in M, and /i a (7£(g)-valued 
semibasic one- form on M satisfying the compatibility condition (23). Let Y be 
the /i-prolongation of order n of A. If F : 5a -^ T5a, we say that A is a 
/i-symmetry of A. 

In the case of scalar equations, the possibility of using /x-symmetries to 
perform symmetry reduction relied ultimately on two facts: (i) the space Xx 
of A-invariant functions is y-invariant for Y a /t-prolongations of A; (ii) the 
standard and the /i-prolongations of A coincide in Ix- This entails that the 
results valid for reduction of an equation A on Tx under standard symmetries 
extend to the case of /^-symmetries. The same holds in the case of PDEs systems. 

Let us first recall that if A = SJ'{d/dx^) + (p^{d/du°') is a vector field in M, 
we denote by Q" := ip'^ — uf^^ its characteristic vector. Then the A-invariant 
manifold in J'^"'>M is the subset Ix C J'^"^M identified by DjQ'' = for aU 
a — 1, ...,q and all multiindices J with < \J\ < n — 1. 
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Theorem 8. In the hypotheses of theorem 6, let Y be the ^-prolongation of the 
vector field X. Then Y coincides with the standard prolongation of the same 
vector field X on 2x ■ 

Proof. This follows from tlicorcm 6 and the definition oiXx- Indeed, write 
^j in the form "^"j = \l/j + Ff, see theorem 6, and suppose that for \J\ = k 
the difference term Ff is written as a combination of the DjQ^, i.e. Fj = 



(r-^)JJDjQ^ Then from (22) we have 



F^, ^ 6^,[D,ir't]iDjQ^) + iA,r,[iT''tiDjQ^)+DjQ'^] ; 

this is again a combination of terms of the form DjQ^. Thus if the Fj vanish 
on Tx for \J\ = k, the Fj with | J| > fc also vanish on Tx- 

Note that F°- = {Ai)'^Q^, so that the condition is satisfied for \J\ ~ 1, and 
the proof of the theorem follows by the recursive computation above. <C> 



6 Examples 

In the examples below we will consider PDEs in two independent variables, 
{x, t). In this case we will also write X = £,dx + rdt + 'fdu, and ji = adx + /3di. 

6.1 yu-symmetries of given equations 

In order to determine /i-symmetries of a given PDE A of order n, one can 
proceed in the same way as for ordinary symmetries. That is, consider a generic 
vector field X acting in M, and its /i-prolongation Y of order n for a generic 
/x = Aidx% acting in J^^^M. One then apphes Y to A, and restricts the obtained 
expression to the solution manifold Sa C J^'^'^M. The equation A, resulting by 
requiring this is zero is the determining equation for /x-symmetries of A; this is 
an equation for ^, t, ip and A^, and as such is nonlinear. 

If we require Xi are a function on J^^' M , all the dependencies on uj with 
\J\> k will be explicit, and one obtains a system of determining equations. This 
system (or the equation A*) should be complemented with the compatibility 
conditions between the A^. 

If we determine apriori the form /i, we are left with a system of linear equa- 
tions for ^, T, ip] similarly, if we fix a vector field X and try to find the ^ 
for which it is a /x-symmetry of the given equation A, we have a system of 
quasilinear equations for the A^. 

The heat equation 

Let us first consider the heat equation 

we will use the ansatz /i = A^dx* (here x^ — x^x"^ — i)-, 

A, = D,P{x,y,u) ; (24) 
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this guarantees that the compatibihty condition (10) is satisfied everywhere (not 
just on Sa)- 

Proceeding as mentioned above, we obtain the determining equations for 
/x-symmetries of the heat equation [under the ansatz (24)]; these resuh to be 



2P„T + 2t„ == , 






2P^T + 2r^ = , 






Pu^ + PuuS. + 2P„e« + en« = , 






P^T + PuuT + 2PuTu + T„„ = , 






-(fixx + Vy- 2(pxPx - (fiPx - ^P'xx + ^Py = , 




PxT + PxxT - PyT + 2PxTx + taUxx -Ty+ 2Px^ + 2£,x = 


-0, 


2PuPxT + 2PxuT + 2PxTu + 2PuTx 


+ 2txu + 2Pu£. + 2^„ 


==0 , 


-ipuu - 2ipuPu - ^Pl - VPuu + 2PuPx^ + 2Pxu^ + 2P. 


^u 1 ^-'tfSa; 



^sa:i/ ^ U 7 



-2(^,„ - 2ipxPu - 2ipuPx - 2ipP^Px - 2ipPxu + -Pic + Pxx^ - Py^- 
1 ^-'a;S2: i ^xx sy '-' ■ 



After some (lengthy but completely standard) computations, we obtain that the 
more general solution to these is given by 



-p 



T{x,t,u) = e^^ 
ip{x,t,u) — e^^ 



ci + C2t + {cz/2)x + {ci/2)xt\ , 
"c5 + cat + (c4/2)i2] , 

C{x,t) + (-(c2/2)x- (c4/8)(x2-2a;t)+C6)■ 



where Ci are arbitrary constants, and C,{x, t) is an arbitrary function satisfying 
Ct = (.xx- Thus, we just obtain the standard symmetries of the heat equation 
[4, 10, 13], with the factor exp[— P(a;,y,u)]; this is no accident, but follows 
from the ansatz (24), see remark 5 (see also [2]). The characteristic Q := 
'P ~ C^a; ~ TUt will be the same as for standard symmetries (with a factor e^^), 
and the symmetry reduced equations will give nothing new. 

The Euler equation 

Let us consider the Euler equation 

ut + uux ^ Q ; 

we will write as usual X = S^dx + rdt + ^du, and /i = adx + f3dy. 

The condition for X to be a y^-symmetry for the Euler equation is that 



ipUx + aUip + f3(f> + u'^UxCtT + UUx/St — UUxOtS, — Uxl3^+ 
+ ipt + UUxTt - Uxit + Ufx + U^UxTx - UUx£,x = . 



2., . „,„, c - n (25) 



This should be complemented with the requirement that DxP — Dta when 
Ut + uux — 0. With the ansatz 

a — a{x, t,u) , P — (3{x, t, u) , (26) 

(note (24) cannot be verified in this case if a and /3 do actually depend on u) 
the dependence of the equations above in Ux is explicit, and (25) splits into two 
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equations: 

{au + f3)ip + ipt + uipx = ; 

(p + {au^ + I3u)t — {au + /3)^ + uTt — Ct + u^t^ — u^x — . 

These are again nonlinear equations for the functions (a, (3,^,T,ip). A special 
solution is provided e.g. by 

a = u , (3 ~ — u^/2 ; 

^ = , T = [B{u) - A{u)t/u] exp[-(uV2)t] , ^ = M^) exp[-{u^/2)t] . 

Note that for this /i, the compatibility conditions Dta ~ D^P is satisfied only 
on the solution manifold S'a, see remark 4. 

This /x-symmetry corresponds to a nonlocal ordinary symmetry Z of expo- 
nential type, see remark 5. We have in facts 

6.2 Equations with given ^-symmetries 

We can also consider the opposite question, i.e. given a vector field X and a 
form /i — Aidx* satisfying (10), determine the equations of a given order n which 
admit X as a /i-symmetry with the given jj,. 

To solve this problem, we have to consider the /i-prolongation F of X to 
j('^)M, and solve the characteristic equation for it. In this way we obtain the 
differential invariants for Y , and any equations which is written in terms of these 
will admit X as a (strong) /^-symmetry. 

Example 1 

As a first example, to be dealt with in detail, we will consider /i-prolongations 
of the scaling vector field 

X — xdx + 2i9t + udu ■ 

The invariant coordinates (j/, v) and the parametric coordinate a in M ~ 
{(x,i,u)} can be chosen as ct = x , y = x^ /t , v = u/x; the corresponding 
inverse change of variables is x = cr , t = a'^ /y , u = av. 

It follows easily that in the symmetry-adapted coordinates, X — crda] hence 
the function v — v{a, y) is X-invariant if and only if Uo- == 0, as required by 
the general method (indeed, by the very definition of symmetry-adapted coor- 
dinates) . 

Applying the procedure described in sect. 4, we have 

Ux=v + 2yvy + ava ; ut = -{y^/(7)vy . 
The above can be inverted to give v^- = (l/x)[ux + 2{t/x)ut — u/x] = —Q/x'^, 



V. 



V 



[t'^ /x^]ut. Similarly, at second order we get 

Uxx = 2uo- + 6{y/cr)vy + (7v„„ + Ayv„y + A{y'^/a)vyy 
Uxt = -3{y^/<7^)vy - {y^/a)v^y - 2{y^/a'^)vyy , 
Utt = 2{y^/a^)vy + {y^/a^)vyy . 
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We will now consider the simplest nontrivial choice for fi, i.e. fi = Ada;, with 
A a real constant. 

The second /^-prolongations can be written in the form (20) (with q = I), 
and general explicit expressions for the coefficients are obtained using either 
(12) or (14). With our choice for ^, one gets 

\]/=^ = A (m - xux - 2tut) , ** = -ut , 

\^xx _ _y_^_^ _ 2X{xUxx + 2tUxt) + A^(?i — XUx — 2tut) , 

*"=* =: -2u^t - X{xu^t + 2tutt + ut) , *" = ~3utt ■ 

By the method of characteristics, we obtain with standard computations 
that the invariants of Y in J^^^'M are given by: 

y := (x2/t) , V := (u/x) ; 



Ci 



— xut , C,2 '■= (ujx — 2tut/x — Ux)e^^ 

= xtutt , m := {xut + 2xtutt + x'^Uxt) e'*'^ , 

— {l/x) [(1 — \x){u — xux) + 2Xxtut + x^Uxx + ^xtuxt + ^f^utt] e^^^ 



Thus, any equation A := F [y, w, Ci, ^2, ?7i, 772, ^l — 0, with F an arbitrary 
smooth function of its arguments, admits the vector field X given above as a /x- 
symmetry, with ji — Xdx. Moreover, for {dF/dC,2f' + [OF/ 81^2)^ + {dF/drj^Y ¥" 
0, X is not an ordinary symmetry of A. 

Finally, it is easily seen by restricting the functions given above to Xx that 
the /x-symmetry reduced equation, providing X-invariant solutions to A, is given 

by 

H[y,vXi,m] := F[y,i;,Ci, 0,7/1, 0,0] = 0. 

Let us discuss a completely concrete example. Consider the equation A := 
rii — C,2 — 0. This is written as 

xtutt + (ux + {2t/x)ut — u/x) e'^^ = 

in the original coordinates; in the adapted ones it reads 

y^Vyy + 2y'^Vy + ae^^v^ = . 

The corresponding reduced equation is y'^[yvyy + 2vy] — 0; the general solution 
to this is v{y) — ci + C2/y, where Ci are real constants. Going back to the 
original coordinates, the corresponding solutions are u{x, t) — {cix^ + C2t)/x. 

Example 2. 

Consider the same X as above, with /i — — (l/t)dt. In this case the y-invariant 
functions are spanned by 

Ci^Ux ; C2 = ut/x - {u/x - Ux)/{2t) ; 

rji = XUxx , m = Uxt + XUxx/{2t) , 

773 ^ {xuxx - 'iux)l{^t^) + Uxt/t + utt/x - ut/{xt) + 3M/(4xt^) . 
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Any equation given by F[y, v, Ci, C2, '?i7 ^2, ^73] = will admit X as a (strong) ^- 
symmetry, and if {dF/dCT)^ + {dF/drj2)'^ + (dF/drj^)'^ 7^ 0, X is not an ordinary 
symmetry. 

The restriction of the invariant functions to 2x yields Ci = ^x, C2 = 0, and 
?7i — xuxx, '72 = ^3 = 0; hence the reduced equation will be simply 



H[y, v,Ci,m] := F[y,t.;Ci,0;77i,0,0] = 



Example 3. 



Any equation of the form i^[j/,w, Ci, C2, 771,772,773] = 0, where F is a smooth 
function of its arguments and we have defined 



y ~ \/x^ + lF , f = u , 

Cl = Ux/X , (2 = Mt - {t/x)Ux , 

771 = {y^/x^){xUxx " Ux) , ?72 = {y/x^){xtUxx " X^Uxt ~ tUx) , 

773 = Utt + {t/x^){xtUxx - 2x^Uxt -tUx) , 

admits the rotation vector field 

X = xdt — tdx 

as a /i-symmetry, with /i = — (l/a:;)da;. 

The functions y, v provide invariants in M , and we can select a = arctg(t/a;); 
the inverse change of coordinates is given by x = ycos{a), t — j/sin(cr), u = v. 
The vector field X is then expressed && X ~ da 

The invariant subset Xx is in this case identified by Ut = (t/x)ux, Uxt = 
{xUxx — Ux){t/x'^), utt — [xt^Uxx + {x"^ — t^)ux]/x^. The restriction of the 
invariant functions to Ix yields Ci = Ux/x, (,2 = 0, 771 = [r^ j x^){xuxx — Ux), 
V2 — 0, ?73 ~ Ux/x (note 773 = Ci); hence the reduced equation will be simply 

H[y,v, CuVi] :- Fiy, z;; Ci, 0; 771, 0, Ci] = 0. 

Example 4. 

In previous examples the functions A and i^ in /i = Xdx + i>dt were always 
depending only on x and t; in this last example they will depend on first order 
derivatives of the u. 

Any equation of the form F{y, v; Ci: C2; 771: ^72, 773) = with F a smooth func- 
tion of its arguments, which are 

y — t , V = u/x ; 

Ci=u + log[l - u/{xux)] , C2 = ~iuut)/{x'^Ux) ; 

771 = -e^" [w/(a;^Mi^)] [u'^{ul — Uxx) + xv?^ — (1 + xux)uu1\ , 

772 = e" [u^/{x'^ul.)\ [xUxUxt - {Ux + XUxx)ut] , 

% = [u/{x^ul)\ [xulutt - 2xUxUtUxt + {2Ux + XUxx)Uf] , 



admits the scaling vector field 

X = xdx 
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as a (strong) ^-symmetry, with /i — u^dx + utdt. 

The invariant set Xx is identified by Ux — u/x, u^t ~ ut/x, u^x ~ 0. 
Restriction of the invariant functions to Tx yields Ci = Ci '■— u/x — v'^, (2 = 
^2 := —Ut/x, 771 — r/2 = 0, 773 = 773 := uu/x. Hence the reduced equation wiU 
be simply 

H[v,v,e2A] ■■= F[y,i;;i;,C2°;0,0,%0] = 0. 

6.3 Systems of PDEs 

Finally, we consider systems of PDEs with a given /i-symmetry, and the cor- 
responding reduction. We will denote independent and dependent variables as 
{x, y) and {u, v) respectively. 

Example 1 

Let us consider the scaling vector field 

d ^ d d ^ d 

X ^ x— + 2y— + u— + 2v— 

ox oy ou ov 

and the form /i ~ XI dx with A a real constant; this corresponds to matrices 
Aj given by A(^x) — ^I and A(j^) = 0. 

By applying the vector /i-prolongation formula (21), or using theorem 6 
and (22), we determine the second /^-prolongation Y oi X . We can then solve 
the characteristic equation for the flow of Y in j'^'^^M, and determine a basis of 
K-invariant functions. Such a basis is provided by the following set of functions: 

p ~ x^ jy , w\ :— u/x , W2 '■— v /x^ ; 



Ci 

Cs 

J?! 
'73 

V5 



= yuy/x , (2 := {ux + 2yuy/x - u/x)e^^ , 
= Uy , C4 := {vx/x - 2v/x^ + 2yvy/x^)e^'' 



2 3 

{xvxy + 2yvyy)e^'' , 

\,Jy LO^IJ I JjUjij I ^ Jj y LL '^(1 11 J C- 



(27) 



Vxx + 4:v/x - 3vx/y + 4:yuxy - '^yvy/x + Ayuy/x+ 
-Ay'^Vyy/x'^ + Sy'^Uyy/x + X{vx - 2v/x + 2yvy/x)\e^^'-' , 
776 := [~Ux + u/x + XUxx + Ayuxy + Ay'^Uyy/x+ 
+\{-u + xux + 2yuy)]e'^^'' . 

Any (system of) second order equation of the form 

F* [y,wi,W2;Ci,---,C4;'7i,---,'76] = (28) 

with F"^ {i — 1, ...,n) a smooth function of its arguments, admits AT as a (strong) 
/Li-symmetry. 

^We stress that this expression for C,\ is not obtained by a direct substitution: indeed now 
(3 = means u^ = u/x\ the general expression for C,i given above becomes singular, but the 
expression for ^^ guarantees that Ux is constant, and actually equal to u/x = d on Q = 0. 
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In order to consider the /z-symmetry reduced equation, it suffices to consider 
the restriction of the functions QiVj on Tx (note that the /i-prolongation and 
the ordinary one coincide on Tx, see lemma 5.) The manifold Tx is identified 
by Q = DxQ — DyQ == 0; in the present case these mean 

Uy ^ {u- xux)/{2y) , u^y = -'{xuxx)/{2y) , Uyy = -{uy + xuxy)/{2y) ; 

Vy = (2w - xvx)/{2y) , v^y = {v^ - xVx^)/{2y) , Vyy = -{xVa;y)/{2y) . 

Substituting these into (27) above, we obtain the expressions for the reduc- 
tion of first and second order y-invariants restricted to Tx , which are 

C'^ = iu-u,x)/{2x), C2° = 0, 

CgO = (2w - xv,)/{2y) , C4° = ; 

ryO = (.t3 {-v^ + xv^^))/{4:y^) , 

jjO ^ (a;3 (^y _|_ a; (u^ _,_ xuxx)))/{'iy^) , 

77O = , r;2 = , r;0 = , r;0 = . 

Thus, the X-invariant solutions of (28) are obtained as solution of the reduced 
system of equations 

W[y,wuW2;C"„Clvi,V2] -^ F^ [y,z«i,z«2;Ci", 0,(3°, O;/,?,//^", 0,0,0,0] =0. 

Example 2 

Consider next the elementary vector field X = xdx, and the form fi = A(j.)da; + 
A(^y^dy corresponding to matrices 

_ /O -y\ . _]_f~xy {x^-x)y^ 
" ' lo ; ' '^^- y^[ 1 xy 

In this case, a basis for y-invariant functions on J^^' M is provided by 

P^y , wi^u , W2^v ; 

Ci = xux - x'^yv^ , C2 = a;W:r , 

Cs = (l/y)[ywiy - a;^2/Wx - x{xux - a:^2/Va;)] , 

C4 = % - (1/2/^) (a^Ua; - x^yvx) log(a;) ; 

771 = x^Uoja; - 2x'^yvx - 2x^yvxx , ^72 = a;^Wa;a; , 

773 = x[Uxy - x{2Vx + aJl^xx + yVxy)] , 

rji = -{l/y'^){xy{xvx + x^Wxx - yvx?/) + x[ux + x{uxx - iyvx - 2xyVxx)] log(a;)} , 
775 = {l/y'^)[Ax'^Ux + 2a:^Mxa; + y{-2x^Uxy - -ix^Vx - 2x'^Vxx+ 

+yuyy - 2x^yvxy + 2x^yVxy)] , 
% = (l/2/^){-2x2y7;x - 2x^yVxx + y^Vyy+ 

+2x[ux + y{xvx + x'^Vxx + xyvxy - Uxy)] log(a:)} . 

Any system of second order equation of the form 

F'[y,wi,w2;Ci,--X4;Vi,---,V6] = , 
_F* a smooth function of its arguments, admits AT as a (strong) /i-symmetry. 
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The system identifying Tx is now given by 

XUx —- u, Ux I Xllxx ~" ^1 -^^xy — ^i 

xvx = 0, Vx+ xvxx = 0, xvxy = ; 

Tx is the linear space on which u^ = u^x — u^y = v^ = v^x = Vxy = 0. 
Restriction of the invariant functions given above to this space yields 

Cl = C2 — 1 C3 — "a J Ca = "^y \ 

rli^ril^ril^ril:^Q , rf^ ^ Uyy , ryg = Vyy . 
The reduced system is therefore 

W[y,U,V;Uy,Vy,Uyy,Vyy] = F"" [ij , U , V ', , , Uy , Vy ; , , , , Uyy , Vyy] = 0. 
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